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ABSTRACT 

We use large-scale simulations to investigate the morphology of reionization during the final, overlap phase. 
Our method uses an efficient three-dimensional smoothing technique which takes into account the finite mean 
free path due to absorption systems, A a b s , by only smoothing over scales R s < A a b s . The large dynamic range of 
our calculations is necessary to resolve the neutral patches left at the end of reionization within a representative 
volume; we find that simulation volumes exceeding several hundred Mpc on a side are necessary in order to 
properly model reionization when the neutral fraction is ~ 0.01-0.3. Our results indicate a strong dependence 
of percolation morphology on a large and uncertain region of model parameter space. The single most important 
parameter is the mean free path to absorption systems, which serve as opaque barriers to ionizing radiation. 
If these absorption systems were as abundant as some realistic estimates indicate, the spatial structure of the 
overlap phase is considerably more complex than previously predicted. In view of the lack of constraints on the 
mean free path at the highest redshifts, current theories that do not include absorption by Lyman-limit systems, 
and in particular three-dimensional simulations, may underestimate the abundance of neutral clouds at the end 
of reionization. This affects predictions for the 21 cm signal associated with reionization, interpretation of 
absorption features in quasar spectra at z ~ 5-6, the connection between reionization and the local universe, 
and constraints on the patchiness and duration of reionization from temperature fluctuations measured in the 
cosmic microwave background arising from the kinetic Sunyaev-Zel'dovich effect. 

Subject headings: cosmology: theory — dark ages, reionization, first stars — intergalactic medium 



1. INTRODUCTION 

On scales relevant to the large-scale structure of reioniza- 
tion - hundreds of Mpc - numerical calculations cannot self- 
consistently track all the sources and sinks of ionizing radia- 
tion. Resolving the IGM on the small scales relevant to QSO 
absorption systems optically-thick at the Lyman edge, also 
known as Lyman limit systems (e.g., Weyman et al. 1981; 
Tytler 1982; Mo & Miralda-Escude 1996), is particularly dif- 
ficult. These absorption systems were likely to have impeded 
the progress of reionization (e.g., Miralda-Escude et al. 2000; 
Gnedin & Fan 2006). 

What are likely values for the comoving absorption system 
mean free path at the end of reionization, A a b s ? Prochaska 
et al. (2009) determined the abundance of Lyman-limit ab- 
sorption systems, dN/dz = 1.9[(1 +z)/4.7] 5 ' 2 , over the red- 
shift range 3.6 < z < 4.4, which translates into a comov- 
ing mean free path of A a b s = c/H(z)/(dN/dz) ^415 Mpc 
at z = 3.7. This dependence on redshift is steeper than re- 
ported previously (e.g., Storrie-Lombardi et al. 1994), and im- 
plies A abs oc (1+z) -6 ' 7 . However, Songaila & Cowie (2010) 
extended the constraints out to z ~ 6, finding a much shal- 
lower dependence, A a b s oc (1+z)" 3 - 44 , with a measurement 
of A a b s ~ 34 Mpc at z ~ 5.7. Comparison of the absorption 
spectra and hydrodynamic simulations can also constrain A a b s , 
with Bolton & Haehnelt (2007b) reporting 20 < A abs < 50 
Mpc at z ~ 6. The evolution of the mean free path at z < 6 
leaves open a wide range of possible values during the reion- 
ization epoch, from several to hundreds of Mpc. Our goal is 
to determine the sensitivity of the history and morphology of 
reionization to this wide range of possible mean free paths. 
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Direct simulation of all the processes involved in reioniza- 
tion is not currently possible. Theoretical studies focus ei- 
ther on the small scale hydrodynamics and radiative transfer 
(e.g., Miralda-Escude et al. 2000; Gnedin 2000; Haiman et al. 
2001; Shapiro et al. 2004; Ciardi et al. 2006; Pawlik et al. 
2009; Raicevic & Theuns 201 1; McQuinn et al. 201 1), or on 
the large scale morphology of reionization (e.g., Abel & Wan- 
delt 2002; Furlanetto et al. 2004; Iliev et al. 2006; McQuinn 
et al. 2007; Trac & Cen 2007; Shin et al. 2008; Friedrich et al. 
2011). On the largest scales, it has been somewhat surpris- 
ing how well simplified "semi-numerical" approaches (e.g., 
Zahn et al. 2007; Mesinger & Furlanetto 2007; Thomas et al. 
2009; Choudhury et al. 2009; Alvarez et al. 2009; Santos et al. 
2010; Zahn et al. 2011a; Mesinger et al. 2011a) match much 
more computationally expensive radiative transfer techniques 
(e.g., Abel et al. 1999; Gnedin & Abel 2001; Sokasian et al. 
2001; Abel & Wandelt 2002; Iliev et al. 2006; McQuinn et al. 
2007; Trac & Cen 2007). In order to efficiently survey the 
wide parameter space of A a b s and focus on the large scale mor- 
phology and overall progress of reionization, we have chosen 
the simplified semi-numerical approach. Our implementation 
is similar to that of Zahn et al. (2007), with the the addition 
of a treatment of photon consumption in absorption systems. 
Our approach, in which we treat abundance of absorbers as an 
input parameter, is complementary to that taken by Crociani 
et al. (2011), where the semi-numerical approach was used 
to determine the large scale distribution of absorbers, rather 
than to model the effect of the absorbers on the progress of 
reionization, as we do here. 

Our results are parametrized by different values of for the 
minimum source halo mass, ionizing efficiency of collapsed 
matter, and the absorption system mean free path. We de- 
scribe our simulations in §2, followed by the main results in 
§3, ending with a discussion in §4. All calculations were done 
using (n m ,fl A ,h, cr 8 , n,) = (0.27, 0.73 ,0.72, 0.8, 0.96), consis- 
tent with WMAP 7-year data (Komatsu et al. 201 1). All dis- 
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TABLE 1 
Simulation parameters 
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"Obtained using sharp i-space filtering 

tances are comoving. 

2. SIMULATIONS 

Table 1 shows our simulation parameters: £, the ioniz- 
ing efficiency, A a b s , the comoving absorption system mean 
free path, f ev , the photoevaporation time for the evolving 
A a bs models, and M m ; n , the minimum halo mass of galaxies. 
Also shown are the redshifts when the global ionized fraction 
equals 0.5, 0.9, and 1 - zo.5, zo.9, and z ov , respectively, and the 
bubble mean free path when the neutral fraction is 0.1, Ao = 
\( x m = 0.1). The parameters were varied to give a Thomson 
scattering optical depth of(r_, to, t + )=(0.06, 0.09, 0. 12), corre- 
sponding to the 2-cr constraint from WMAP (Komatsu et al. 
2011), with the He I fraction tracking H I, and instantaneous 
He II reionization at z = 3. All simulations used 4096 3 cells in 
a 2 Gpc//i box. 

2.1. Mean reionization history 

We model the mean reionization history wtih two spatially 
uniform mean free paths: that corresponding to ionized bub- 
bles, Ab, and Lyman-limit systems, A a b s - Ionizing radiation 
is attenuated by their superposition, so that A~' fp = A^ 1 + A~b S . 
The spatially averaged ionizing flux is F(z) = A m f p (z)e(z), 
where e(z) is the ionizing photon emissivity. We set e(z) = 
C n H,o/coii(zX where £ is number of ionizing photons per col- 
lapsed hydrogen atom, corrected for recombinations outside 
of absorption systems. The mean ionization rate is 



dt °° A abs (z) + A b (x)' 



(1) 



When A a bs ^ Ab, photons typically reach the edges of bub- 
bles without being absorbed by intervening Lyman-limit sys- 
tems, and the photoionization rate is equal to the emission 
rate. When A a bs <C Ab, the probability of a photon reaching 
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FIG. 1. — Bubble mean free path, Ab, as a function of mean neutral frac- 
tion, .vhi. All models exhibit a characteristic scale of ~ 10 Mpc at the half- 
ionized epoch. During percolation, lower values of A a b s lead to lower val- 
ues of Ab. The solid line corresponds to the evaporating minihalo model 
with f C v = 100 Myr. At jchi < 0. 1 the relation approaches approximately 



Ab oc x. 
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the edge of a bubble is A a bs/Ab, and the ionization rate is sup- 
pressed. 

A more realistic treatment would involve allowing F, A/,, 
and Aabs to vary spatially, obtaining a spatially-dependent so- 
lution of equation (1), which can then be averaged to find the 
global reionization history. We have chosen the simpler uni- 
form model of equation (1) as a starting point. Although this 
choice does not affect our results on the morphology of the 
ionization field at fixed ionized fraction, the assumptions un- 
derlying equation (1) should be kept in mind when interpret- 
ing the reionization history and photon consumption rates we 
find. 

2.2. Spatial variations 

Our three-dimensional model is based on that of Furlan- 
etto et al. (2004), later extended to three dimensions by Zahn 
et al. (2007). Its main assumption is that a region is fully ion- 
ized if its collapsed fraction is greater than some threshold, 
C/coii > 1. As shown in Alvarez et al. (2009), by smoothing 
the linear density field over a range of scales, one can effi- 
ciently determine when each point is first reionized, z r - We 
do not smooth over scales with radii larger than A a bs, since 
absorption systems shield radiation from these distances. We 
assume that the absorption systems contribute a spatially uni- 
form opacity, and therefore the scale beyond which we do not 
smooth is the same everywhere. 

This approach is desirable because the order in which points 
are ionized is the same as that in which the Furlanetto et al. 
(2004) criterion is first met around each point. However, the 
reionization redshifts do not generally result in the correct av- 
erage ionization fraction for a sharp smoothing filter in real 
space (see below for solutions using a sharp filter in A:-space, 
which does conserve photons; Zahn et al. 2007). To obtain 
a self-consistent solution, we determine a correction to the 
reionization redshift at each point, z r Az r ), by matching the 
simulated reionization history to the global reionization his- 
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FIG. 2. — Global reionization histories (bottom) and photon consumption (top). Different colors indicate different values of A a b s , while the panels indicate 
different Thomson scattering optical depths, t cs = 0.06, 0.09, and 0.12, from left to right. Note that in all cases the early part of the reionization history is 
insensitive to A a t, s , when the bubbles are relatively smaller. Lower values of A a j, s delay the end of reionization. 



tory, 



dp 
dz 



dz : 



dx dt 
dt dz 



dz, 



(2) 



where dp/dz is the distribution of uncorrected simulated 
reionization redshifts, and dx/dt is given by equation (1). 
This calibration procedure does not change the overall topol- 
ogy we find at fixed ionized fraction, but does change the 
overall reionization history. 

We calculate X\,(x) by random raytracing. If a ray starts in 
an ionized region, it is extended in a random direction until 
it reaches a neutral cell. The bubble mean free path is then 
the average over the lengths of all such rays. The average 
converges with on the order of 10 4 rays. 

About 30 logarithmically-spaced smoothing scales are suf- 
ficient to converge on the reionization morphology. Because 
we determine the reionization redshifts simultaneously, the 
amount of data stored and the operation count are greatly re- 
duced compared to determining the ionization field at each 
redshift. Since the assumption that points are either highly- 
ionized or neutral is made in both cases, no information is 
lost by following the more efficient approach. 

2.3. Sharp k-space filtering 

Using a sharp A:-space filter leads to a method which does 
conserve photons (Zahn et al. 2007) in the long mean free 
path limit, at the expense of using an unphysical filter which 
"rings" in real space, and an ambiguity in the assignment of 
a smoothing scale Rf to the filter kf - particularly troubling 



given that we want to exclude scales Rf > A a b s - However, for 
comparison we have calculated a few cases with sharp &-space 
filtering, shown in Table 1 . 

The ionization history is found by integrating a diffusion 
equation for the distribution in density of points that have not 
crossed the barrier (Bond et al. 1991), 



dIL(A,S,z) 
dA 



1 <9 2 n d5 on 

2 ~df~dA ~dS' 



(3) 



with boundary condition LT(A,0,z) = 0, where S = B(A,z)- 
5, 5 is the density contrast, B(A,z) = 6 c (z) — erF (1 — 

C" 1 )^^ -A)is the bubble barrier (Furlanetto et al. 2004), 



and cr^ in is the variance of the z = linear density field at the 
scale of M m i„, with the growth of structure encoded in the time 
dependence of 5 c (z). 

We smooth only over scales larger than A a b s by starting in- 
tegration of equation 3 at A,- = (T 2 (A a b s ) - the variance on the 
scale of the mean free path is matched to the one obtained 
using a top-hat filter in real space - with initial condition 



n(A,-,M = 



exp 



(^-g(A,-, z )) 2 

2A; 



(4) 



The total ionized fraction is obtained by integrating over the 
distribution at A = er 2 • , 



x(z)=l- 



n(cdj.,tf,z>/5. 



(5) 
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If the initial condition for II is specified at A = 0, (i.e. starting 
at a scale — > oo), then this corresponds to smoothing over all 
scales, and 



x(z) = Cerfc 



8 e Q0/yJ: 



2al 



■■ C/ce>n(z)- 



(6) 



2.4. Evolving A a bs due to minihalos 



To model an evolving value of A a b s , we assume that the 
absorption systems are minihalos(e.g., Abel & Mo 1998). 
Given that the halo cross section increases as M 2 / 3 , while 
dn/dlnM oc M~ l , small objects should dominate the mean 
free path, subject to photoevaporation effects (e.g., Haiman 
et al. 2001 ; Shapiro et al. 2004). 

Consider a uniform distribution of dark matter halos with 
a number density n/,(z) = dn/dlnM/,. We will assume that in 
neutral regions, halos at mass M/, have gas fractions of unity 
and would serve as Lyman-limit absorptions systems, while 
in ionized regions, such halos survive for a time f ev before 
being evaporated by the ionizing background. The number of 
absorption systems evolves according to 



dn 



;ibs 



dz 



dlnx 

dz 



(«h-«abs) + £ev(l+z) 5/2 «ab S (z) 



(7) 



where £ ev = (ffo^m^ev)" 1 ^ 260 (* ev /100 Myr)- 1 . We choose 
M h = 5.7 x 10 3 M Q [(1 +z)/T0] 3//2 , corresponding to the cos- 
mological Jeans mass, and consider three models, in which 
f ev = 10, 50, and 100 Myr. These timescales are relatively 
long compared with the rather short photoevaporation times 
for low-mass halos found by Shapiro et al. (2004). However, 
if reionization is "photon-starved", as observations are sug- 
gesting (e.g., Bolton & Haehnelt 2007b), then the flux could 
be low at the end of reionization, leading to longer evapora- 
tion times. 
The mean free path is 



(8) 



where R vlI corresponds to the halo virial radius. However, the 
boundaries of the absorption systems are not well-defined due 
to hydrodynamic effects and departures from spherical sym- 
metry. For our purposes, the minihalo model we describe here 
is sufficient to incorporate the effect of an evolving absorp- 
tion system mean free path in our calculations. More accurate 
modeling of absorption systems during reionization will need 
to incorporate the radiative transfer of ionizing radiation in 
cosmological hydrodynamics simulations which resolve the 
Jeans mass. 

We start with an initial guess for A a b s (z), only allowing 
points to cross the barrier at some scale and redshift if that 
scale is below the current value of A a b s . The values obtained 
for Xb(x) from the simulated ionization field are then fed back 
into equations 1, 7, and 8, and the process is reapeated until 
convergence is reached. 

3. RESULTS 

Fig. 1 shows the bubble mean free path as a function of 
the neutral fraction for several parameter choices. At the 
half-ionized epoch, all models exhibit a characteristic bub- 
ble scale of Ab — 10 Mpc, with models with rarer sources (i.e. 
increasing efficiency of ionizing radiation) leading to some- 
what larger bubbles. During percolation, most of the varia- 
tion in Ab(xei) comes from variation of A a b s , which is most 
pronounced at jchi < 0.2. The large values of A;, indicate that 
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FIG . 3 . — (top) Neutral fraction vs. redshift for the three evolving A a b s mod- 
els considered here, in addition to the reference case, for which \ bi — > oo 
(right-most dotted curve), (bottom) Absorption mean free path vs. redshift. 
As the last neutral patches are cleared away, so are the last reservoirs of 
cold neutral gas in which the minihalos can survive as absorption systems, so 
A a b s increases at the end of reionization. As the evaporation time decreases, 
the destruction rate of absorption systems increases, leading to longer mean 
free paths. Also shown are the constraints on A a b s from Bolton & Haehnelt 
(2007b) at z = 5 and 6. 

simulation volumes in excess of several hundred Mpc on a 
side are necessary in order to properly model the last ten per 

cent of the reionization process. At x m < 0.1, A& ~ x m > in 
agreement with the model of Miralda-Escude et al. (2000), in 
which a constant number of neutral clouds decrease in size at 
the same fractional rate at the end of reionization. 

The reionization history for all of the fixed A a bs models is 
shown in Fig. 2. The early evolution is not sensitive to A a b s , 
because the photon mean free path is determined by the size 
of the ionized bubbles themselves. At later times, when the 
bubble sizes become comparable to A a b s , the reionization his- 
tory is sensitive to A a b s , as can be seen by comparing the evo- 
lution for the bracketing values of 8 Mpc/h and 256 Mpc/h. 
For the A a b s = 8 Mpc/h case, the redshift at which the overlap 
occurs is delayed by Az ~ 1.5, relative to A a b s = 256 Mpc/h, 
for r es = 0.09 (see also the "z ov " entries in Table 1). 

Also shown in Fig. 2 is the number of ionizing photons per 
ionized atom, 



» 7 fe) _ C/coll(z) 
«HIl(z) X(z) 



= 1 + 



1 



z} AbQQ 
x(z) Jo A abs (x') 



dx' . 



(9) 



As Aabs — y oo, x — y C/coii and n 7 /«Hii ^ 1- We only count 
photons up to the moment of overlap (x = 1), as we are con- 
cerned here with the consumption of ionizing photons dur- 
ing the reionization process. This integral converges, since 
at the end of reionization we expect Xb(x) oc (1 -x)~ 2 ^ as the 
remaining diffuse neutral hydrogen patches dissappear. For 
A a bs = 8 Mpc/h, about 3 photons per atom are consumed in 
the absorption systems by z ov . In the longest mean free path 
case, A a bs = 256 Mpc/h, the corresponding fraction is about a 
half. 

Shown in Fig. 3 is the evolution in neutral fraction and ab- 
sorption system mean free path for the three minihalo ab- 
sorption system models that we simulated. Models with 
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longer evaporation times lead to a higher abundance of ab- 
sorption systems, and hence shorter mean free paths and a 
relative delay in the time of overlap. The mean free path in 
the ? ev = 100 Myr model evolves from about 40 Mpc when 
xm ~ 0.1, to about 80 Mpc at overlap. This model is also 
plotted as the solid black line in Fig. 1, which shows that the 
instantaneous Ab values for the evolving A a bs model roughly 
interpolate between those for the fixed A a bs models. This can 
also be seen by noting that A () = A b (x H i = 0.1) for the evolving 
mean free path model, for which A a b s (x H i = 0.1) ~ 32 Mpc/h, 
is the same as that for the corresponding fixed mean free path 
model with A abs = 32 Mpc/h, A = 170 Mpc. 

Fig. 4 shows the ionization field in a 5-Mpc//i slice at 
x = 0.75 for M min = 10 8 M Q , with A abs = 256 Mpc/h (top) and 
A a bs = 8 Mpc/h (bottom). The two cases are quite different, 
with many more small neutral patches for A a b S = 8 Mpc/h rel- 
ative to A a bs = 256 Mpc/h. Also shown in Fig. 4 is the reioniz- 
tion redshift in a Q.S-Mpc/h slice. Lower values of A a b s result 
in a more extended reionization overlap period - i.e. the max- 
ima in reionization redshift are higher, and the minima are 
lower. 

4. DISCUSSION 

We have carried out large-scale simulations of reionization 
in a 2 Gpc/h volume, including a finite mean free path to 
absorption systems. Absorption systems have a significant 
effect on the characteristic scales at the end of reionization. 
For the Thomson scattering optical depth reported by WMAP, 
T es ~ 0.09, we find that the characteristic bubble size when 
the universe is 90 per cent ionized is quite sensitive to A a b s , 
with A — 440 Mpc for A abs = 256 Mpc/h and M min = 1O 9 M , 
while A ~ 94 Mpc for A abs = 8 Mpc/h and M^ D = 1O 8 M . 

Calculations using a sharp fe-space filter lead to a more 
modest extension of the percolation phase (see Table 1) com- 
pared to the reionization history obtained by the solution of 
equation (1), on which the reionization histories shown in Fig- 
ure 1 are based. For r es ~ 0.09, the difference in the redshift 
when the universe is 90 per cent ionized between the A abs = 
8 Mpc/h and A abs = 256 Mpc/h cases is only Az .9 ~ 0.3, 
while the corresponding case using a sharp real space filter 
and equation (1) is Az .9 ~ 0.7. The delay in the overlap time 
in both cases is Az ov ~ 1.5. 

Using either sharp fe-space smoothing or equation (1) to 
model the effect of absorption systems on the reionization 
history has drawbacks. While equation (1) conserves pho- 
tons in the long mean free path limit, a mean flux and opac- 
ity are assumed, and further work will be necessary to vali- 
date and/or improve upon the approach layed out in section 
2.1. Using a sharp fc-space filter results in a reionization his- 
tory which also conserves photons in the long mean free path 
limit (k F — ► 0), but at the expense of using an oscillatory fil- 
ter which can "leak" photons from high density regions into 
lower density ones, and for which there is no unique choice 
for the relation between A a b s and kf . While our results for the 
timing and photon consumption of the end of reionizaiton are 
obtained from equations (1) and (9), it is important to note 
that more work will be necessary to test their accuracy. How- 
ever, results for the morphology of reionization (Figures 1 and 
4) do not depend on the model for the global history, and are 
therefore more robust. 

Our results are consistent with those of Furlanetto & Oh 
(2005), in which consumption of ionizing photons in dense 
systems extends the end of reionization considerably. Choud- 
hury et al. (2009) used the semi-numerical approach in a vol- 



ume 100 Mpc/h across, and found a similar trend with de- 
creasing A abs , indicating a transition to an "outside-in" mor- 
phology near the end of reionization. We find that neutral 
patches may also remain in voids, where formation of ioniz- 
ing sources is delayed and radiation from the nearest sources 
is shielded. Previous studies that modeled the physical origin 
of absorption systems as minihalos (e.g., Ciardi et al. 2006; 
McQuinn et al. 2007) also found similar effects to those we 
find here, although these studies were more focused on the 
intermediate stages of reionization and on smaller simulated 
volumes, where the photon mean free path was not as small 
relative to the bubble sizes as we find in our low mean free 
path cases in the final, percolation phase. 

As mentioned in the introduction, our approach is comple- 
mentary to that taken by Crociani et al. (2011), in which the 
semi-numerical approach was used to determine the distribu- 
tion of absorbers during reionization. They found that their 
spatial distribution is quite inhomogenous, owing both to in- 
trinsic density fluctuations in the IGM, as well as fluctuations 
in the ionizing radition field - regions far from sources have 
a relatively low flux, resulting in a higher abundance of ab- 
sorbers. This affect was also pointed out by McQuinn et al. 
(2011), who used high-resolution hydrodynamic simulations 
post-processed with radiative transfer to determine the mean 
free path as a function of flux, finding a strongly nonlinear de- 
pendence of the mean free path on the emmissivity of ionizing 
radiation. 

These effects imply an important improvement to our 
model is not only a time-varying mean free path, but also a 
spatially varying one. To accomplish this, more work will 
need to be done on the "sub-grid" physics during reionization, 
using the results of high-resolution cosmological simulations 
with radiative transfer of a background radiation field, cou- 
pled to the hydrodynamics of the gas. However, given that 
the final patches to be ionized are those most distant from the 
most luminous sources, the delay in the very end of reion- 
ization that we find here is likely to persist when the inho- 
mogeneity of the IGM opacity due to absorption systems is 
properly taken into account. 

The lingering neutral clouds we find would further com- 
plicate interpretation of quasar absorption spectra at z ~ 6 
(e.g., Bolton & Haehnelt 2007a; Lidz et al. 2007; Alvarez 
& Abel 2007; Wyithe et al. 2008; Furlanetto & Mesinger 
2009). As discussed by Mesinger (2009), gaps in quasar spec- 
tra can come from either these neutral clouds or mostly ion- 
ized but optically-thick absorption systems. The line-of-sight 
abundance of these neutral clouds at fixed ionized fraction in- 
creases with decreasing absorption system mean free path. 

As shown by Weinmann et al. (2007) and Alvarez et al. 
(2009), the Local Group may have been reionized by external 
sources, i.e. the progenitors of the Virgo Cluster at a distance 
of ~ 20 Mpc. If A a b s < 20 Mpc, however, local reionization 
would have been delayed until Local Group prongenitors had 
formed in sufficient number (e.g., Munoz et al. 2009). Since 
the satellite population is sensitive to the timing of reioniza- 
tion (Busha et al. 2010), lower values of A a b s could imply a 
higher local satellite abundance. 

Finally, our results have important implications for the in- 
terpretation of temperature fluctuations in the cosmic mi- 
crowave background at multipoles I ~ 3000 (Das et al. 
2011; Reichardt et al. 2011), through the kinetic Sunyaev- 
Zel'dovich (kSZ) effect (e.g., McQuinn et al. 2005; Iliev et al. 
2007). In particular, Mesinger et al. (2011b) used nearly the 
same three parameters and approach as we used here for their 
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FIG. 4. — left: Reionization redshift in a 0.5 Mpc//i-thick slice for the A a bs = 256 Mpc//i simulation (top) and the A a b s = 8 Mpc//i simulation (bottom). Both 
cases were for t cs = 0.09 and M m i n = 10 8 Mq. Clockwise beginning from top-left, regions shown are 2 Gpc/h, 500 Mpc/A, and 125 Mpc//i across. The radii of 
the circles correspond to A a |j S . right: Same as left panels but for the projected ionized fraction in a 5 Mpc//i-thick slice at a time when the mean volume-weighted 
ionized fraction is 0.75. Note that the ionization field contains considerably more small scale structure in the A a b s = 8 Mpc//i case. 



parameter study of the dependence of the kSZ signal on the 
reionization scenario. Reducing the mean free path shifts 
power to higher multipoles because of the accompanying de- 
crease in the bubble size that we find. Although the kSZ ef- 
fect is only just beginning to constrain the duration of reion- 
ization (Zahn et al. 201 lb), more detailed future observations 
and analysis will begin to provide constraints on the patchi- 
ness in addition, and understanding the effect of absorption 
systems on reionization will be crucial. 
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Toronto. This work was partially supported by NASA ATFP 
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